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TEMPERATURE DISTRIBUTION OF A THIN-WALLED, 
TRANSPARENT, SPHERICAL, EARTH SATELLITE 


by 

Joseph T. Skladany 
Goddard Space Flight Center 


INTRODUCTION 

Investigations of the temperature distribution of spherical satellites began almost twenty years 
ago. As early as 1955, Sandorff and Prigge (Reference !) wrote a paper that included conduction in 
the skin but excluded radiation within the hollow sphere. The sphere was considered as heated by 
solar input only, earth input being neglected. The skin was divided into a number of “stations” or 
“nodes,” and the temperature distribution was obtained by solution of a number of simultaneous 
equations. A paper on the general thermal design of satellites by Schach and Kidwell (Reference 2) 
made use of a similar subdivision method; this nodal technique is still in use today. Hrycak (Refer- 
ence 3) considered both skin conduction and internal radiation effects and concluded that, for thin- 
walled satellites, internal radiation plays the dominant role. 

Wood and Carter (Reference 4) gave equations for the mean temperature and hot- and cold-spot 
temperatures on a thin-walled spherical shell at the sub-solar position and at a position in the earth’s 
shade. Conduction was neglected, and an approximate method was used to account for albedo inputs. 
In the computation of the local body temperatures, distribution of the earth inputs over the skin of 
the satellite was ignored. It is generally acknowledged that, for a sphere at 1000-km altitude, earth- 
emitted and albedo energy incident on part of the sphere can be as much as three times the average 
over the entire sphere; this may drastically affect local temperatures. 

A proposal (Reference 5) to NASA suggested the use of a partially transparent sphere of poly- 
propylene with an etched pattern of aluminum on both sides to increase reflectivity. The temperature 
at a few special points on the sphere and the mean temperature were obtained from tracings of the 
infinite number of reflections at the inner surface. It was this proposal that prompted the present 
investigation. 

This paper is concerned with temperature distribution over the surface of a thin-walled, trans- 
parent sphere in an earth orbit. The only mode of heat transfer considered is radiation, since lateral 
and radial conduction in the skin can be shown to be negligible (Reference 3). It is assumed that the 
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thermal mass is small and that transient effects are therefore negligible. The sphere is assumed to be 
a gray body; only specular transmittance is considered. This paper is not concerned with thermal- 
property measurements; it is assumed that the materials have the required properties. 

Before making a thermal analysis of a transparent sphere, we examine the case of an opaque 
sphere in order to introduce the theory and display the necessary thermal inputs. Once the expression 
for the temperature distribution of an opaque sphere is derived, attention is focused on a transparent 
enclosure. The heat transfer equation for radiation, including an infinite number of internal reflec- 
tions, is solved by Poljak’s radiosity method, described by Eckert and Drake (Reference 6), This yields 
three simultaneous integral equations which, by means of the thermal model technique (Reference 2), 
can be solved for any shape enclosure. Spherical geometry is then introduced, and a solution for the 
mean radiation temperature of a sphere is obtained. Finally, an expression for the temperature at any 
point on the sphere is derived. Both expressions are found to be general in that they yield mean tem- 
perature and temperature distribution for any sphere, whether opaque or transparent, or for a sphere 
having opaque portions. The solution of this last case can be derived from the effective thermal prop- 
erties (Reference 5). 

THE OPAQUE SPHERE-MEAN TEMPERATURE 

General Discussion 

The mean temperature of a satellite in equilibrium with its environment is found when the ther- 
mal energy absorbed by the satellite is equated to the thermal energy radiated by the surface of the 
satellite* 

Hie external thermal inputs of an earth satellite consist of radiant energy from the sun, solar 
radiation reflected by the earth (albedo), and radiation emitted by the earth (Figure 1). Since the solar 
radiation and albedo extend over the entire solar spectrum, the fraction of this energy absorbed by the 
surface is given by its absorptivity. Also, since the earth-emitted energy lies almost wholly in the infra- 
red portion of the spectrum (as from a 250 K blackbody. Reference 7), the fraction of this energy 
absorbed by the surface is given by its emissivity. 

The rate of solar energy impinging on an element normal to the sun's rays is the well-known solar 
constant (0.0333 cal/sec-cm 2 ) and at any other orientation is simply a function of the cosine between 
the solar vector and the normal to the surface (Reference 2). If this input is integrated over the sur- 
face of a sphere, the total solar energy impinging on the sphere equals the solar constant multiplied by 
the projected area of the sphere. The albedo- and earth-radiation inputs to a sphere can likewise be 
found as functions of the projected area of the sphere. The energy radiated by the satellite, however, 
is a function of its entire surface area. The energy balance on the sphere can be written in integrated 
form as in Reference 2, 


+ % P o A P + e o ^o A p Aoe Q 



(I) 
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FROM SUN 



Figure i -Energy incident on an earth satellite. 


where 



(For explanation of symbols, see Appendix C— Nomenclature.) This equation defines the satellite 
mean temperature for radiation exchanges, or, solving Equation 1 for T m , we obtain 



If the satellite is far removed from the earth, this reduces to 



Earth Inputs 

In the previous section, earth radiation and albedo were introduced, but no mention was made as 
to how they are determined. In this section, earth radiation will be discussed. Albedo will be dis- 
cussed in the following section. 
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To at least a first approximation, the earth can be regarded as a blackbody radiating diffusely at 
250 K (Reference 7). The earth-radiated energy directly incident on the satellite’s surface can there- 
fore be derived from the equation 


G 0 = oT e F l- 2> < 4 > 

where I and 2 refer to the satellite and the earth, respectively, as shown in Figure 2 and F^ 2 the 
shape factor from area 1 to area 2. To compute Fj_ 2 accurately it would be necessary to set up an 
expression for the shape factor from a point on the satellite’s surface, integrate over that portion of 
the earth viewed from the point, and finally integrate over the satellite’s surface. However, if two sur- 
faces obey Lambert’s law, the shape factor from a point on one surface to the whole of the other is 
merely the solid angle subtended by the second surface at the point on the first surface. Since the 
earth is much larger than the satellite and the distance between them is great in comparison to the 
satellite’s diameter, the solid angle subtended by the earth at any point on the satellite disk is very 
nearly a constant. Hence, it is sufficient to consider the satellite a disk and compute one shape factor 
for a point on the satellite disk (as shown in Figure 2) and multiply this factor by the projected area of 
the sphere to obtain the fraction of the earth-radiated energy that impinges on the sphere. 

By definition under the above assumption, 

_ _ iff cos 0 cos <pdidk 

f <u-a d - F 1-2 ~ J J y > (5) 

p A d \ 


where 

/ = distance between the satellite center and an element on the earth’s disk, 

0 = angle between / and the normal to an element on the earth’s disk (or satellite disk), 

A d - area of earth’s disk visible from the satellite, 

and 

A p = projected area of satellite. 

However, di cos 0// 2 is the solid angle dco subtended by the satellite disk at a point on the earth’s disk. 
When the satellite’s diameter is small compared to the distance between the earth and satellite, dai is 
very nearly a constant for any point on the earth’s disk and is equal to A p /l 2 . With this substitution, 
Equation 5 becomes 


F 


1-2 



cos (pd k 


( 6 ) 


By reference to Figure 2, we obtain 


r cos 2 6 
cos 0 


4 




Figure 2— Earth radiation to a sphere. 


and 


dk = 2 wl sin sin 0}] = 2m 2 cos 4 0 tan 0 sec 2 0of0 , 

where 0 is the angle between the line from the satellite's center to the earth's center and the tangent 
from the satellite's center to the earth's surface* Substituting these expressions into Equation 6 and 
simplifying, we obtain 




- F 


1“2 



2 sin 0d0 = 2(1 - cos 6) , 


( 7 ) 


but 


cos 8 - 



( 8 ) 


When we define h = rfR e7 he,, when h is in the dimensionless form of earth radii. Equation 8 becomes 


^ 1 - 2=2 


1 - 


(A 2 - 1) 


1/2 


( 9 ) 
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Finally, 


G 0 = 2oT* 


1 - 


ih 2 - 1) 


1/2 


(A plot of this expression is shown in Figure 8 in a later portion of this report.) 


( 10 ) 


Albedo 

The earth-reflected solar radiation incident on a spherical satellite can be obtained in much the 
same manner as the earth radiation. The reflectance of the earth varies over a wide range, with the 
mean ranging from 0.34 to 0.5 (Reference 8), depending on the portion of sky covered by clouds. 

The problem is set up as in Figure 2, but now the flux from the earth is the reflected solar energy 
instead of the earth’s radiated energy, and the integration must be limited to the sunlit portion of the 
earth; this complicates the limits of integration. This problem has been investigated in some detail by 
Cunningham (Reference 9); we shall make use of his results. To do this we must know the altitude of 
the satellite and the angle between the earth-sun line and earth-satellite line. The altitude is known 
from orbital predictions; the angle can be computed from latitude and longitude observations of the 
sun and satellite, as is explained in Appendix A. 


THE OPAQUE SPHERE-TEMPERATURE DISTRIBUTION 
General Discussion 

The temperature distribution of an opaque spherical satellite can be found by means of an energy 
balance set up at a differential element on the sphere’s surface. Again, the inputs to the external sur- 
face of the sphere consist of direct solar energy, albedo, and earth radiation. The net energy leaving 
the internal surface of the element must also be considered, since this is not necessarily the same at all 
internal surface elements. A schematic of the thermal energy balance is shown in Figure 3. The 
energy balance equation at any element di is 

a 0r^0i + a oi^oi + e oi G oi " oe oi + ' (H) 

where dQ t is the net, inward-directed energy leaving the internally facing surface of di. Using Poljak’s 
radiosity method (Reference 6), we obtain 


dQ ( = (B. - H.)di , 

(12) 

H n 

+ + 
& ^ 

(13) 
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Figure 3 —Thermal balance on an opaque element 
of an earth satellite. 


and 


and, substituting 13 in 14, 


diH. 


■I 


B dk-di B k dk 



B dk~df e k + Pp^P Bk ’ 


( 14 ) 


05 ) 


(Differential element di is not unique; an equation written for di is true for all other elements.) 

Since at any instant the sphere is at equilibrium with its environment, the total energy incident 
on the inside of the sphere is equal to the total energy leaving the internal surface of the sphere; that 
is, in equation form. 



dQ ( 




H k dk . 


After B k is substituted from Equation 13, Equation 16 becomes 


or, rearranged. 


I. 

i 



(e k p*H k )dk , 


0 -PlW k dk. 


( 16 ) 


07 ) 


( 18 ) 
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Next we define F dk _ di , the shape factor from element di to element dk. 

The geometry for this problem is shown in Figure 4. By reciprocity, F dk-<u dk = F di-dk dL fi y 
definition. 


F = 
r di~dk 


cos cos <j> k dk 

id 2 


( 19 ) 


where 

4>. = the angle between the normal to area i and the line joining area i to area k, 
4> k = the angle between the normal to area k and the line joining area k to area i, 
and, from geometry, 

, , 1 / 
cos <pj — cos 

where r' is the radius of the satellite; therefore, 

_ fidk = dk 
t di-dk = 4ir 2/ 2 l 2 A ' 

Substituting this expression in Equation 15, we obtain 

»,-j/ <«*+<>;»»><»• 


( 20 ) 


( 21 ) 
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If p% = constant p *, Equation 18 can be combined with Equation 21 to give 


**r 



e k dk 

1 -p* ' 


Substituting for H. in Equations 12 and 13 yields 


dQ i = e ( di - di 



( 22 ) 


( 23 ) 


but 



= eoT* 

m 


and finally 


dQ. = oe(T*-T*)di; (24) 

Equation 1 1 becomes 

Vo,- + a Q P 0j + e Q G 0( . = a(e 0 . + e t )T* - oeT * . (25) 

Therefore, to find the temperature at any point on the surface of the sphere we need only define 
the external inputs at that point and solve Equation 2 for T . The next section deals with defining 
earth inputs. 

Earth Inputs to Differential Element on Sphere 

To define the earth inputs to a differential element on the satellite, it is necessary to find a shape 
factor from the element to the earth’s surface (sunlit portion of the earth for the albedo) (refer to 
Figure 5). Again, by definition, 


] f cos a 1 cos a. 
f dhdk nj ~ D ~ 


dk 


( 26 ) 


where 

= the angle between the normal Nj to the element on the satellite and the line joining the 
element to an element on the earth’s surface, 


a 2 = the angle between the normal N 2 to an element on the earth and the line joining this ele- 
ment to the element on the satellite’s surface, 

D - the distance between the element on the satellite and the element on the earth. 
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The integration takes place over that part of the surface of the earth which is visible from di (and also 
illuminated by the sun in the case of albedo). 


From vector analysis we obtain 


cos O-j 


Nj ■ D N 2 ■ D 

is^iDi md cos “2 = is^5r- 


(27) 


These expressions can be written in terms of 6, <p, and These, when substituted in Equation 
26, yield an unwieldy integral that requires numerical integration. Rather than attempt this, we will 
make an approximation that allows use of existing data. 


Approximation for Earth Inputs 

The approximation for earth inputs consists of imagining that each element of the earth is located 
at the sphere’s center, as shown in Figure 6. Because of the relative sizes of the earth and satellite, 
locating an element in this manner has very little effect on D t a y , or a 2 as long as the element’s orienta- 
tion remains unchanged. It is now permissible to use existing literature on earth-radiation and albedo 
inputs to a small flat plate (References 7 and 10). 

To do this we must find X, the angle between the surface normal and the earth-satellite line, and 
(for the albedo) 6 s , the angle between the earth-satellite line and the earth-sun line. (Expressions for 
these angles can be found in Appendices B and A, respectively.) 



Figure 5— Earth inputs to an element 
on a sphere. 


Figure 6— Approximation for earth inputs 
to element on a sphere. 
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THE TRANSPARENT SPHERE 


General Discussion 

In finding the temperature distribution of a transparent sphere, we deal with external inputs in 
the same manner as with an opaque sphere. However, with a transparent sphere, part of the energy 
passes through the surface and impinges on the inner surface. The problem of accounting for this addi- 
tional energy will be discussed in this part of the report. 

Besides the assumptions made previously, we will also assume that energy is neither scattered nor 
diffracted as it passes through the surface of the sphere, This means that the path of a ray of energy is 
not changed as it passes through the skin of the satellite. We find the external energy that impinges on 
the interior skin of the satellite by multiplying the external energy that impinges on the exterior skin 
by the appropriate transmissivity. 

Poljak’s radiosity method (Reference 6) is again extensively used in this section and enlarged to 
include solar energy and albedo as well as infrared energy. 

Tire thermal balance is treated generally at first, so that the method may be used for any shape 
enclosure. Late'r, the spherical shape is introduced. This allows some simplifications to be made. 


Net Internal Energy Leaving an Element of an Enclosure 


Figure 7 shows the thermal balance on a transparent element of a satellite. Note that this figure 
includes the following energy terms: /, the solar energy incident on the interior skin; L, the solar 
energy leaving the interior skin; r/, the solar energy transmitted through the element; and t*H, the 
infrared energy transmitted through the element. By redefining dQ. so that it includes not only net 
energy leaving the interior skin but also transmitted energy, we can see that Equation 1 1 still holds for 
a transparent surface. We can now write 


Total energy incident 1 _ + j . 

on an interior element f ‘ 


(28) 



Figure 7 -Thermal balance on a transparent element 
of an earth satellite. 
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Since some of this energy passes through the surface, we can write 

Total energy leaving j 

the element (includ- >= di(B i + L. + T j I i + r*7/ f ) , (29) 

ing transmitted) ) 

where the last two terms take into account that incident energy which passes through the surface. 


Subtracting Equation 29 from Equation 28 gives 


dQ i - di[B i +L r { 1 - r f )/ f - (1 - rf)//.] . 

(30) 

From the definitions of the terms in Equation 30, we obtain 


B i = e i +P*H., 

(31) 

L i = Pi I i> 

(32) 

diI i = di(S j +P i ) + J dkF^L,, 

(33) 

and 

r 


diM; = diG,. + / dkF dk . m B k . 

i t 

(34) 

Substituting Equations 3 1 and 32 in Equation 30 gives 


dQ t = di\[e t + <p* + rf - 1 )//,.] + (p- +r ( - l)/.j 

(35) 

and, using reciprocity, dkF dk _ dj = dlF dhd r Substituting Equations 3 1 and 32 
respectively, gives 

in Equations 34 and 33, 

/,-«S ( +J> ( )+ / F^'l 

(36) 

and 

r- 


H, = G,+ (e k +f>tH k )F d ,. dk . 

J k 

(37) 


Equations 35, 36, and 37 give three simultaneous equations at a differential element. These equa- 
tions are completely general and therefore hold for an enclosure of any shape. Since the integrals are 
difficult to solve analytically, we can most easily solve these expressions by subdividing the surface 
into a number of finite nodes, changing the integral signs to summations, and solving the resulting set 
of equations simultaneously. This network method is similar to the method discussed by Schach and 
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Kidwell (Reference 2). It is usually true that accuracy increases with the number of nodes chosen; 
obviously solution very soon becomes unwieidly and requires a computer. 

When the surface in question is a sphere, some simplifications can be made, resulting in a closed- 
form solution without the need to subdivide the surface into a number of finite nodes. This will be 
investigated in the remainder of the report. 

Mean Temperature 

To find the mean temperature of a transparent sphere, integrate Equation 1 i over the sphere’s 
surface. Then the left-hand side becomes the same as in Equation 1 ; the right-hand side must now in- 
clude the net energy leaving the interior of the sphere. When dQ. is integrated over the surface we 
obtain 


Q = A[e + (p* +t* - \)H + (p + r - 1 )/] , (38) 

where the properties are constant over the surface. 

From the definitions of H and /, from the assumption that the transmitted energy passes through 
the surface with no diffraction, and from the fact that the shape factor from the interior of the sphere 
to itself equals unity (since the sphere is a closed boundary), Equations 36 and 37 become 

rA pi s o + 'V 

Al=TA p (S 0 +P 0 ) + l>AI= (39) 

and 

AH = T*A p G Q +AB . (40) 

If we integrate Equation 30 over the surface of the sphere with p* constant, B i becomes B\ we can 
substitute for B to obtain 




1 -p* 


(4!) 


Hence the energy balance becomes 


(p* + r* - l)(r*A G q + Ae) 

a 0^p^0 + + e O A p G 0 “ A ° € Q^m + Ae + l -o* 


(j> + t-1)tA p (S 0 +P q ) 
+ 1 "P 


(42) 


Noting that 


and 


P* + t*- l=-a* = -e 

p + t- 1 = -a 
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and rearranging terms, we rewrite Equation 42 as 


+ v(“o + pf^)* V^'o + 7~;) . (43) 

which reduces to Equation 2 for an opaque sphere. If the inner surface has the same properties as the 
outer surface, Le., a = a Q and e = e then 

"W + v (‘ + r^) + e ' 4 o c o( l + iy) = <*< (‘ + r^) ■ (43a) 


Temperature Distribution 


As on the opaque sphere, the temperature at any point on the transparent sphere can be obtained 
from Equation 11. Note that dQ. must be derived from Equations 35, 36, and 37. 

With constant surface properties, and making use of the expression for the shape factor between 
two elements on a sphere (Equation 20), we can rewrite Equations 36 and 37 as 


and 


/. = 0S ( .+ J P.) + p/ 

(44) 

//. = c. + e + p*H . 

(45) 


If substitutions from Equations 39 and 41 are made for / and //, Equations 44 and 45 become 

/. = S + P- + p 

i i ( H 


V (5 0 + 'o)| 

44(1 -p) j 


and 


H i -G i + £? + p 4 


T*(A!A)G 0 +e 


1 ~P' 


Therefore Equation 35 can be rewritten as 

dQ, 


di 


= e - e 


e + 




prA( s +P) 


and, simplified. 


IF ^,-ocs, 

- <T *@’) (t^ c <> ■ ■ eT ( 1 + F^) e • 


(46) 


(47) 


(48) 


(49) 
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EARTH RADIATION * 10*^ 


By the use of Equation 43, one of the input terms in Equation 49 can be eliminated. Since <7 fl 
is usually the least important term, we will choose to retain this term. This allows the user to simplify 
the result, when G Q becomes negligible (at altitudes above approximately 1000 km), by neglecting this 
term (see Figure 8). 


Solving Equation 43 for ^(5^ +P Q ) and substituting in Equation 49 yields 


d Qj_ 

di 


e r a(S j +l\)-€(G j )-~f(G 0 ) 


epV 

1 - 0 *' 



- oTi 



(50) 



Figure S— Earth radiation incident per square centimeter projected area of a sphere vs altitude. 
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Finally, substituting Equation 50 in Equation 1 1 and solving for Tf yields 


T t ~ a(e '+ e ){ a o (S o > + p o ii> + a ^ S i + P i ) + e o G o, + eG / 


&V e Jhl apr ( e ° i-p*) 

W^ 0 l-P* „ J ar\ 

1° 1-pj 


U-P)| 


e + e„ 




1-P* 


(51) 


EXAMPLES 

To illustrate the expressions derived, a fictitious orbit will be used* The temperature distribution 
in the skin of a satellite will be found for the case of a satellite at the subsolar position in orbit (ie., 
where the sun, satellite, and earth are in line) and also 180 degrees from this orbital position, where 
the only input is earth radiation. These positions were chosen because they are both points of sym- 
metry (Le., distribution of temperature is only in the polar direction) and because the first case gives 
the hottest condition and the second gives the coldest, both cases being important in thermal design. 

For our fictitious orbit we will assume that the satellite is launched into a circular, polar orbit 
along the prime meridian at 12 noon GMT on the first day of spring. This positions the sun on the 
equator and in the plane of the orbit. From Appendix A, for these conditions, the sun’s position is 
given by 


S= i 

and, when at the subsolar point, the satellite’s position is given by 

L = i ; 


therefore 



From Appendix B, 


n = i sin 6 cos 0 + j sin 8 sin 0 + k cos 9 ; 
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therefore 


cos P = sin 6 cos 0 
cos X = - sin 0 cos 0 . 


r 


and 


However* since we have located the satellite at the subsolar point we can set 6-90 degrees and 
vary 0 from 0 to i8G degrees to cover all points on the sphere* because of symmetry. Then 

cos = cos 0 


and 


cos X = - cos 0 , 


or 

and 


P-0 

X = 1 80^ “ 0 . 


Knowing d s and X* we can now refer to the literature (References 7 and !0) and find the albedo 
and earth radiation to every point on the exterior of the satellite. Also, since the unit internal normal 
is _ n and because of our assumption of no diffraction, the interior incident albedo and earth radiation 
can be found in the same manner if we let X = 0 and multiply the input by the appropriate transmis- 
sivity. 

Also note that the solar input can be obtained from Sdi cos P for 0° 90° and that the solar 

input equals zero for 90° < p < 1 80° * For the satellite at the subsolar point, this becomes Sdi cos 0 
for 0° < 0 < 90° and zero for 90° < 0 < 180°. Likewise, the solar input to the interior of the sphere 
is - rSdi cos ( 1 80° - 0) for 90° < 0 < 1 80° and zero for 0 < 0 < 90°, 


The albedo and earth radiation incident on the exterior of the sphere are plotted in Figure 9 as a 
function of 0. We can use the same Figure for internal inputs, substituting 0 for (180° - 0) and multi- 
plying albedo values by r and earth-radiation values by t*. 


The albedo on the entire sphere at an altitude of 1000 km and 6^ = 0 from Cunningham (Refer- 
ence 9) is 44,8 W/cm 2 X 10^ 3 or 10,7 X IQ -3 cal/sec-cm 2 , and the earth radiation from Figure 8 is 
5,28 X 10" 3 cal/sec-cm 2 . 


At 1 80 degrees from the subsolar point the only input is earth radiation, which is assumed con- 
stant over the orbit (5.28 X 10“ 3 cal/sec-cm 2 ). 


To fully illustrate the equations derived, three spherical satellites will be used, two opaque and 
one transparent. The surface properties used are summarized in Table 1. 

The temperature distributions are shown plotted as TjT m versus 0 in Figures 10 and 1 1 . 

A sample computation for the polypropylene sphere is shown at both positions in orbit. The 
mean temperature and temperature distribution are computed, and the inputs and results for 0 < 0 < 
180 are tabulated in Table 2 for the subsolar point and in Table 3 for shade position. 
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Figure 9 “Albedo and earth-radiation distribution incident on a spherical satellite 
at the sub solar point and 1000- km altitude. 



Table 1— Thermal properties. 


Sphere 

a o 

a 

e o 

e 

P 

P* 

T 

T* 

Source of 
measurement 

Black 

1.0 

1.0 

1.0 

1.0 

0 

0 

0 

0 

„ 

Alodine 

0.324 

— 

0.185 

0.70 

- 

— 

0 

0 

GSFC 

Polypropylene 

0.042 

0.042 

0.103 

0.103 

0.05 

0.036 

0.908 

0.861 

Reference 5 




Figure 1 0 “Temperature distribution 
at subsolar position. 


Figure 1 1 -Temperature distribution in earth’s shade 
at ISO degrees from subsolar position. 


In earth shade, (*S 0 +P 0 ) = 0; therefore 

T 4 K5.28X10’ 3 ) 

m Aa ~ 4 ( 1 .3 6 )( 1 0~ 12 > 


= 9.7 X 10 8 

and 

T m = 176 K. 


19 




Table 2- Inputs and results for polypropylene sphere at subsolar position. 



Table 3— Inputs and results for polypropylene sphere in earth-shade position. 
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At the subsolar point and <P = 0°, by Equation 51, 


T t = 7"){ a o (5 o< + V + + + e 0 G 0i + *<?, 




ep 


v 4 e ° + r^.) 


t a 

m 

L 

(e + £T * \ 

l 0 l - p7 

e + e 0 

|(l-p>(a„+i-p) 


1 -p* 


or, when numerical values are substituted, 


"4 - 


1 X 3 0' 


- 3 


* (1.36)(1Q -12 )(0. 103 + 0. 103) 


0.042(33.3 + 0) + 0,042(0 + 7.65) + 0. 103(0 + 3.42) 


+ ^(5.28)(10* 4 ) 


(42.0X1 0 8 ) 


/ (0.103X0.861) 

(0.042X0.05X0.908) (0, 103 + — — 

(0.103)(0036)(0.861) \ (1 -0.036) 




(1 - 0.036) 


(3 


0.05) (o. 


042 + 


(0042X0908) \ 
(1-0.05) ) J 


(0.103 + 0 103) 
Tf = 77.1 X 10 8 , 


(0.042X0.05X0908) 0.103 + 


(o. 


(0 103X0.861) 


! ) 


(1 


- 0.05) ^0. 


042 + 


(1 -0036) / (0103)^ 

’ + 1 - 0.036 


(0.042)(0.908) 
(1-0.05) 


! ) 


© 


- = 1.84, 


and 


T i 

-jr= S-165 

1 m 
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In the earth-shade position for <j> = 0 , all the S and P terms are zero, and hence, 


j\4 = 

' <1.36)(10“ 12 X0.103 + 0.103) 


v(0.103)(39.7)(10" 4 ) 


+“(5.28)(10 -4 ) 


9.7(1 0 8 ) 


{0.103)(0.036)(0.861) 


(1 - 0.036) 


/ (0. 103)(0,861)V 

(0.04 2)(0.05)( 0.908) fo. 103 + ~ 

5 . 


(1-0.05) 0.042 + 


(o, 


(0.042X0.908)' 
(1 -0.05) 


(0.103 


(0.042X0.05)10.103 + 


(° 


(0.103X0.861) 


5 


(i- 0,036) / (0.1 03) 2 


+ al03, L (l-0.05)(a04^ ( -5^M«08)\ (1 - 0.036) 


+ (1-0.05) ) 


or 


r 4 = 15.0 x io 8 . 

Hence, 


\V 


= 1.55, 


or 


1.116 
' m 


We treat the opaque cases in the same manner, noting that P., G jt r, and r* are all equal to 
zero, and thereby reduce Equations 43 and 51 to the simpler forms (Equations 2 and 25) derived for 
the opaque sphere. 


CONCLUDING REMARKS 

A method for determining the temperature distribution of a transparent or partially transparent 
nonconducting satellite has been derived. Equations 11,35, 36, and 37 can be solved by numerical 
methods for any satellite configuration with any combination of transparent or opaque surfaces. The 
only restrictions are that the external thermal inputs and the surface thermal radiation properties of 
the satellite must be known. 

By means of the derived expressions and their application to a sphere with constant surface prop- 
erties, expressions for mean radiation temperature and temperature distribution were obtained. These 
are closed-form expressions, easily solved by hand calculation, but their verification in the laboratory 
is quite difficult; indeed, no such experiments are mentioned in the literature. 
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The expression for the mean temperature of a transparent sphere is identical with the expression 
derived in Reference 5. The expressions for hot- and cold-spot temperatures, however, do not agree. 
Apparently this is because only average earth inputs were used, as suggested by Wood and Carter (Ref- 
erence 4), rather than inputs distributed over the surface of the sphere as in the present report. The 
use of average earth inputs would lower hot-spot temperatures, especially in orbits of altitude below 
1000 km. Also, it appears that not all the reflections were accounted for in the ray- tracing technique 
described in Reference 5, 

It is of interest to note the differences in the distribution curves between the opaque and trans- 
parent cases. From Figures 10 and 1 1 it can be seen that the position of the sphere 180 degrees from 
the surfs rays is much cooler than the position directly under the sun’s rays for the opaque cases, 
while the temperatures at these two positions are more nearly equal in the transparent case, This raises 
the possibility of using transparent surfaces as a means of achieving more uniform temperatures over 
the spacecraft. However, use of a completely transparent skin docs not completely eliminate temper- 
ature gradients (sometimes it increases them, as in Figure 10), but it shifts the positions of maximum 
gradient (see Figure 11), Possibly by some combination of opaque and transparent surfaces, nearly 
uniform temperatures may be realized. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland, February 26, 1970 
861-51-75-01*51 
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Appendix A 


Angle Between Earth-Satellite Line 
and Earth-Sun Line 


It is assumed that the satellite’s longitude and latitude are known for any time t. In the following 
equations (see Figure A-l) north latitude and east longitude correspond to positive angles, while south 
latitude and west longitude correspond to negative angles. The American Ephemeris and Nautical 
Almanac gives the latitude (or declination) of the sun for any date. 

The coordinates are set up as follows: the i vector, in the equatorial plane, passes through the 
Greenwich meridian; the k vector, perpendicular to the equatorial plane, passes through the North 
Pole; the j vector, in the equatorial plane, is perpendicular to i and k. Thus a right-hand system is 
formed. 

The direction of the satellite and sun can now be indicated by unit vectors from the earth’s 
center. 



Figure A-l— Geometry of earth-satellite and earth-sun line. 
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At any time t, given in Greenwich Mean Time on a given day of the year, the following informa- 
tion is available: 

For the satellite, 

latitude and longitude, from predictions or sighting. 

For the sun, 

declination, from Reference 1 1 for date given; 

latitude, from time at which satellite is at given latitude and longitude. 

Then, referring to Figure A-l, we obtain 

L = cos (lat) cos (long)L + cos (lat) sin (long)} + sin (lat)k (A-l ) 

and 

S= cos (decl)fcos (/- 12)1 5] i + cos (decl)fsin (t- 12)15]) + sin (decl)k , (A-2) 

where 

lat = latitude of satellite at time t in degrees, 

long = longitude of satellite at time t in degrees, 

decl - sun’s latitude on the given date in degrees, 
and 

/ = Greenwich Mean Time in hours (0 to 24). 

Then, 

L • S 

cos 8 = „ = cos (lat) cos (long) cos (decl) [ cos ( 1 2 - /) 1 5 ] 

s |L||5| 

+ cos (lat) sin (long) co$ (decl) [sin ( 1 2 - r)l 5] + sin (lat) sin (decl) . (A-3) 

It is important to note the sign convention mentioned earlier. 
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Appendix B 

Satellite Coordinate System and Relations 
Between Surface-Normal and Solar Vector and 
Between Unit Normal and Earth-Satellite Line 


In general, most satellites are spin-stabilized; hence the spin axis would be the natural choice for 
one of the coordinate axes. It would then be necessary to transform this coordinate system to a geo- 
centric coordinate system to find the angles necessary for determining earth input and albedo. How- 
ever, since we are considering a non-spinning satellite, we are free to use any coordinate system we 
choose and hence will use a geocentric system, thus eliminating the need for any transformation. 
Figure B-l illustrates this. For certain applications (such as earth-oriented satellites) this selection of 
coordinates may prove unsatisfactory. In such cases the reader is referred to Swalley’s report (Refer- 
ence 12) which deals with this problem in much greater detail. 

With satellite-centered coordinates, the unit normal at any point is given by 


_ W 
" |VF| 


(B-l) 


where F(x,y, z) = 0 is the equation of the surface. 

Since the surface under consideration is spherical, 

F{x,y,z') = x 2 + y 2 + z 2 + R 2 = 0, (B-2) 

where R is the radius of the sphere. Change to spherical coordinates, as shown in Figure B-2 ; now VF 
becomes the vector R and |VF| = R, which gives the unit normal 

n = i sin & cos <j> + j sin 8 sin <j> + k cos d , (B-3) 


where 

0 = angle between k and the normal to the element 
and 

<f> = angle measured in the ijth plane from i to the element (see Figure B-2). 

Then 0, the angle between the normal to the surface and the solar vector, can be computed from 

cos 0 = n • S ; (B-4) 

also X, the angle between the normal to the surface and - L, the satellite-earth line, can be computed 
from 

cos X= n • (-L) = -n ■ L . (B-5) 
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Figure B-I -Geocentric coordinate system. 



At certain positions in the orbit, the satellite may be shaded by the earth. For the satellite to be com- 
pletely in the earth’s shade, two conditions must be satisfied: the satellite’s height above the earth’s 
center multiplied by sin 6 S must be less than the earth’s radius R e , and values of cos 6 S must be nega- 
tive. These conditions can easily be checked at any satellite position. 
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Appendix C 


Nomenclature 


A = surface area of satellite 

A D = area of disk on the earth subtended by the satellite 

A p = projected area of satellite 

B = infrared energy leaving interior surface of satellite per unit surface area 

Bj(B k ) - infrared energy leaving interior surface of area i(k) per unit surface area 

D = vector from an element on the earth’s surface to an element on the satellite’s surface 

- shape factor from area i to area k 

F ' dhdk ~ shape factor from element di to element dk 

G q = earth-radiated energy incident on external surface of satellite per unit surface area 

G Qi = earth- radiated energy incident on external surface of area / per unit surface area 

= earth- radiated energy incident on internal surface of area i per unit surface area 

H - infrared energy incident on interior surface of satellite (including reflections) per unit surface 
area 

- infrared energy incident on interior surface of area i{k) (including reflections) per unit sur- 
face area 

I = solar energy incident on interior surface of satellite (including reflections) per unit surface 
area 

/ f (/ A ) = solar energy incident on interior surface of area i(k ) (including reflections) per unit surface 
area 

L = solar energy leaving interior surface of satellite per unit surface area 

L.(L k ) = solar energy leaving interior surface of area i(k) per unit surface area 

L = unit vector from earth to satellite 
Nj = vector normal to element of satellite 

N 2 = vector normal to element of earth 

P ; - earth-reflected solar energy (albedo) incident on internal surface of area i per unit surface 
area 

P 0 - earth-reflected solar energy (albedo) incident on external surface of satellite per unit surface 
area 

= earth-reflected solar energy (albedo) incident on external surface of area i per unit surface area 
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R = satellite radius vector 

R = radius of earth 

5\ - direct solar energy incident on the internal surface of area / per unit surface area 

£ 0 = direct solar energy incident on the external surface of the satellite per unit surface area 

S 0i = direct solar energy incident on the external surface of area i per unit surface area 

S = unit vector from earth to sun 

T - absolute temperature of earth 

7 \ = absolute temperature of area i 

T m = root mean fourth absolute temperature 
di = differential surface element of area i 

dk = differential surface element of area k 

dQ i - net, inward-directed energy leaving the internally facing surface of di 
e “ emissive power of internal satellite surface = ueT A 

e i = emissive power of internal satellite area / = oe £ T f 

h - r/R e (see below) 

/ = distance between two elemental areas 

n » unit vector nonnal to element of satellite 

r - distance from satellite’s center to earth’s center 

r = radius of satellite 

a = mean solar absorptivity of in ternal surface of satellite 
a. = solar absorptivity of internal surface of area / 
a Q = mean solar absorptivity of external surface of satellite 
a 0/ = solar absorptivity of external surface of area / 

j3 = angle between normal to elemental area on satellite and solar vector 

e = mean infrared cmissivjty of internal surface of satellite 

€. = infrared emissivity of internal surface of area / 

Cq - mean infrared emissivity of external surface of satellite 
e Qi ~ infrared emissivity of external surface of area i 
& s - angle between earth-sun line and earth-satellite line 

X - angle between satellite surface normal and earth-satellite line 

p - mean solar reflectivity of infernal surface of satellite 
p i - solar reflectivity of internal surface of area i 
p* = mean infrared reflectivity of internal surface of satellite 
p * = infrared reflectivity of internal surface of area / 
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a = Stefan- Boltzmann constant 

r = mean solar transmissivity of satellite 
t * = solar transmissivity of area / 

T* " mean infrared transmissivity of satellite 

t * = infrared transmissivity of area / 

cj = solid angle 
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